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INTRODUCTION1 

The  fixed  charge  problem  was  formulated  by  G.  B.  Dantrig  and  W.  Hirsch  In  1954  [l]. 
It  arises  in  situations  which  involve  the  plant, ing  of  several  Interdependent  activities, 
some  or  all  of  which  have  set-up  charges  (or  fixed  charges  independent  of  the  activity 
level  as  long  as  it  is  positive)  associated  with  them.  The  problem  may  be  formulated 
as  follows: 

Min  f(x)  -  0 (x)  +  Z(x) 
subject  to  Ax  ■  b 

*>°  0) 

where 

n 

ZM  -  l  c.x. 

J-l 


and 


0(x)  -  Y  <*j('  - 

j-t 


) 


1  if  *,  ■  0 
J 


0 


if  Xj  >  0 


Am  x  n  1  bm  x  I  ’  Cl  x  n  "  (cl . cn}’  dl  x  n  “  (dl . dn}  are  9’ven  real  matrices  a"d 

x^  i*R°  •  Corresponding  to  any  feasible  solution  x  ,  D(x)  is  known  as  the  fixed  charge 

component  of  the  cost  and  Z(x)  the  variable  cost. 

G.  B.  Dantzig  and  W.  Hirsch  have  shown  that  min  g(x)  is  attained  at  an  extreme 
point  of  the  convex  polyhedral  set  determined  by  (l)  (See  [1,2]).  [3,4,5]  discuss  some 

approximative  algorithms  for  solving  the  fixed  charge  problem,  especially  when  the  under¬ 
lying  structure  of  the  restrictions  (l)  is  of  the  transportation  type. 

*  The  author  is  indebted  to  Professor  R.  Van  Slyke,  Mr.  P..  Chandrasekaran,  and  Profes¬ 
sor  Alan  S.  Manne  for  their  suggestions  and  criticisms. 
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The  algorithm  described  in  this  paper  applies  in  general  to  any  fixed  charge  problem. 
Since  only  the  extreme  points  of  (l),  which  are  finite  in  number,  have  to  be  scanned,  it 
leads  to  the  optimal  solution  in  a  finite  number  of  steps.  However,  the  algorithm  works 
efficiently  when  (l)  is  nondegenerate  and  the  range  in  the  value  of  Z(x)  for  feasible 
x  is  large  compared  to  the  fixed  charges. 

Algorithm  for  the  Fixed  Charge  Problem  Assuming  that  the  Vertices  of  (l)  can  be  Ranked 
in  Increasing  Order  of  the  Variable  Costs  Z(x) 

An  algorithm  for  ranking  all  the  vertices  of  (1)  in  increasing  order  of  the  linear 

functional  Z(x)  is  given  in  Section  2.  For  an  application  of  this  algorithm  we  have  to 

assume  that  min  Z(x)  is  finite.  Here,  min  indicates  that  the  minimization  is 
Cx|(i)]  [x|(D] 

over  all  x  satisfying  (1). 

Case  1:  min  ?(x)  ■  . 

-  Ex | ( l )  ] 

Assuming  that  d  is  finite,  it  Is  clear  that  Z(x)  is  unbounded  below.  Hence  by  Prob¬ 
lem  19,  page  ll«6  of  [6],  we  know  that  min  f(x)  -  -«  iff  the  system  of  equations 

Ex  |  (i )  D 

Ax  ■  0 

x  >  0  (2) 

cx  <  0 

has  a  feasible  solution. 

Hence  in  all  subsequent  discussions  we  shall  assume  that  min  Z(x)  >  -®  and 

MO)] 

hence  that  min  £(*)  >  • 

Lx  |(i )  J 

Case  2:  min  ?(x)  >-“•=»  min  Z(x)  >  . 

-  MO)]  MO)] 

So  in  this  case  it  is  possible  to  rank  all  the  extreme  points  of  (l)  In  Increasing  order 
of  Z(x). 
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Let  Sr  S2 . be  such  3  ranking  and  let  Z(Sk>  -  Zk  .  Then  we  have 

^1  =  ^2  '  ^Ct  =  ”  *1  anb  *et  ^k  C  ‘  Let  °0  be  a  *ower  bound  on  the 

fixed  charge  component  of  the  total  cost  at  any  vertex  of  (l);  i.e.,  <  Dk  \/k 

method  for  obtaining  Op  is  discussed  at  the  end  of  this  section.  The  efficienty  of 

the  algorithm  improves  with  the  nearness  of  Dn  to  the  greatest  lower  bound  of  D  wk 

0  k 

Suppose  we  have  determined  some  Sr  .  Then  it  is  clear  that  the  optimal  solution  to  the 
fixed  charge  problem  must  be  one  of  the  vertices  S  ^  ....  .  w..ere  kr  is  such  that 


and 


k  +1 
r 


Z  >0  -  D„ 

r  r  0 


(3) 


for  any  k  >  k^  we  have 

Zk  +  Dk  -  Zr  +  (Zk  -  Zf)  +  0k  >  Zr  +  Dr  (0k  -  Dq)  by  (3)  And  since  we  know  by  the 
choice  of  D0  that  Dk  >  0Q  yk 

Z.  +  0.  >  Z  +0  for  k  >  k 
k  k  r  r  r 


Hence  it  is  not  necessary  to  rank  all  the  extreme  points  of  (l)  to  solve  the  fixed  charge 
problem.  As  soon  as  S|  is  found  we  get  an  upper  bound  on  the  extent  of  the  values  of 
Z(x)  to  which  we  may  have  to  carry  on  the  ranking  by  using  the  above  result. 

In  general,  suppose  we  have  determined  Sr  .  Then  it  may  be  necessary  to  rank  the 
extreme  points  of  (l)  to  the  extent  that  Z(x)  <  Zf  +  Df  -  .  Now  the  stages  In  the 

algorithm  can  be  described. 

Sta^e^r:  In  this  stage  Sj . have  been  determined.  Let 

min  rA.  ♦  0.  -  Dft] 
k-l . r  K  K  u 


i 

\ 


* 

r 


If  Ar  =  6r  ’  t^ie  a,9°rithm  terminates  and  the  optima)  solution  is  given  by  the  extreme 
point  corresponding  to 


Min  [Z  +  0  ] 
k-1 . r  K  K 

If  Ar  <  6r  >  then  it  is  necessary  to  determine  Zr+j  and  proceed  to  stage  r+)  . 

However,  we  know  that  the  ranking  algorithm  has  to  be  carried  on  to  find  out  vertices  of 
(l)  only  to  the  extent  of  Z(x)  <  Z|  ♦  6r  •  The  various  stages  of  the  algorithm  may  be 

represented  geometrically  as  follows: 


- - - - ¥ - 

Extent  of 

.  ,_T  |  ...  .... 

values  of 

Z(x)  to 

which 

¥ 

"  ranking 

1 

may  have 

to  be 

1  1 

■  .  -  .  carried 

zM 


l 

I 

I 

I 

• 

I 

I 

O 

o 


4 


I 

I 

• 

I 

I 

I 

o 
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I 

I 
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I 

I 

I 

I 

c 

o 

I 


I 

4 


4 


out . 


2  3 
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The  lines  parallel  to  the  k-axis  on  the  diagram  indicate  the  extent  of  the  value  of  Z(x) 
to  which  ranking  may  have  to  be  carried  out.  At  each  stage,  the  horizontal  line  nearest 
the  k-axis  applies.  This  limit  is  improved  at  each  stage. 


To  Find  0^  ,  a  Lower  Bound  on  the  Fixed  Charge  Component  of  Cost  at  anv  Vertex  of  (1) 
Suppose  the  variables  Xj . are  arranged  in  increasing  order  of  the  value  of 


dj  ,  i.e., 


dl  ^  d2  —  £  dn 


Each  extreme  point  of  (l)  consists  of  m  basic  variables,  but  some  of  them  may  be  at 
zero  values  if  the  problem  is  degenerate.  If  we  know  that  there  is  no  degeneracy  in 
the  problem  we  can  take 


0_  “  d,+. . .+d 
0  1  m 

Even  if  (l)  is  degenerate,  it  may  bt  possible  to  obtain  a  lower  bound,  mf  ,  on  the  number 
of  basic  variables  which  are  positive  at  any  vertex  of  (1).  If  (l)  is  not  totally  dege¬ 
nerate  (i.e.,  b  i*  0) ,  then  none  of  its  canonical  equivalents  can  be  totally  degenerate 
anyway.  This  may  help  us  to  get  some  lower  bo-’nd  for  ni|  . 

If  the  constraints  of  (I)  are  of  the  transportation  type  It  is  very  easy  to  deter¬ 
mine  ni|  .  Then  if  all  dj  >  0  ,  we  can  take 

0  1  m. 


A  crude  value  for  0.  is  0  when  all  d.  >  0 
0  j  - 

Corollary  1:  The  algorithm  works  equally  efficiently  if  we  replace  D(x)  be  any  con¬ 
cave  function  and  0-  by  a  lower  bound  to  min  0 (x)  . 

°  rx|(D] 

Then  f(x)  ■  Z(x)  +  D(x)  is  also  a  concave  function  and  It  is  a  well-known  result 
that  the  minimum  of  a  concave  objective  function  over  a  convex  set  occurs  at  an  extreme 
point . 
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1 1 .  An  Algorithm  for  Ranking  the  Vertices  of  (1)  In  Increasing  Order  of  Z(x) 


Consider  the  linear  programming  problem  in  its  standard  form 


min  Z  ■  cx 


subject  to  Ax  *  b 
x  >  0 


We  shall  assume  that  this  problem  has  a  finite  optimum,  i.e.,  that  min  Z(x)  >  -»  . 

Cx  |  ( 1 ) ] 

Then  it  is  well  known  that  there  exists  a  vertex  of  (1)  which  is  optimal  for  the  above 
problem. 

The  algorithm  developed  here  is  an  extension  of  the  simplex  algorithm.  It  helps  in 
ranking  the  basic  feasible  solutions  of  (I)  in  increasing  order  of  Z,  after  the  optimal 
is  obtained  by  the  simplex  method.  It  uses  only  one  step  pivot  operations. 

Let  the  letters  B  and  T  ,  with  any  subscripts  or  superscripts  if  necessary,  de¬ 
note  basic  feasible  solutions  of  (1).  Suppose  in  any  basic  feasible  solution  B  the 

variables  x  . x  are  basic.  We  shall  indicate  this  by 

1  m 

x  cB  i  *  I , . . .  ,m 
r. 

i 

and 

3  ■  fx  , . . . ,x  } 
r.  r  1 

m 


Here  we  are  defining  a  basic  feasible  solution  by  the  set  of  variables  which  are  basic 
in  it. 

Let  S, ,  S  denote  the  minimal  cost  (w.r.t.  Z (x ) )  basic  feasible  solution  and 
1  max 

the  maximal  cost  basic  feasible  solution  respectively.  We  have  assumed  that  S|  exists. 

Consider  any  basic  feasible  solution  B  .  Corresponding  to  any  nonbasic  variable 
Xj<B  ,  let 


7 


g 

C.  ■  the  relative  cost  coefficient  of  the  nonbasic  variable  x.  cor- 
J  J 

responding  to  the  basic  feasible  solution  B  . 


0 


B 

J 


the  value  with  which  the  nonbasic  variable  Xj  enters  the  basis 
in  the  canonical  form  of  the  basic  feasible  solution  B  . 


the  new  basic  feasible  solution  obtained  by  pivoting  on  the  column 

of  x.  in  the  canonical  form  of  B  . 

J 


From  the  simplex  algorithm 


Z(T?)  =  Z(B)  +  0?  C? 
J  J  J 


The  basic  solutions 
B  . 


T8 

J 


for  j 


such  that 


x.rfB 

J 


are  adjacent  vertices  of  the  vertex 


However,  when  (I)  is  degenerate,  several  basic  feasible  solutions  may  repre¬ 
sent  the  same  vertex  and  all  the  adjacent  vertices  of  this  vertex  are  given  by  the  basic 
B 

feasible  solutions,  Tj  ,  corresponding  to  the  various  canonical  forms  B  which  represent 
the  same  vertex. 

Let  11(B)  denote  the  set  of  all  the  adjacent  vertices  of  B  with  cost  value  not 
less  than  that  of  B,  i.e., 


0(B)  ■  fT8  vj  such  that  x^B  ,  C8  >  0}  . 


w 


It  can  be  seen  that  the  canonical  form  corresponding  to  any  of  the  adjacent  vertices  of 
B  can  be  obtained  by  pivot  operations  on  the  canonical  form  of  B  .  If  (l)  is  nonde¬ 
generate,  then  corresponding  to  each  vertex  of  the  polyhedron  there  exists  a  unique 
basis  B  which  represents  it,  and  equation  (4)  holds  for  each  individual  basis. 

However,  if  B  is  a  degenerate  basic  feasible  solution  of  (l),  let  Vg  denote 
the  vertex  represented  by  it.  Let  Bj,...,Br  be  all  the  basic  feasible  solutions  of 
(I)  that  represent  the  same  vertex  Vg  .  Then  we  should  replace  equation  (k)  by 
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0(vB) 


U  {T.P  yj  such  that  x.ifB  and  C.p  >  0]  , 

p=l  J  J  P  J  “ 


1 


(4a) 


where  0(B)  of  (4)  and  0{Vg)  of  (4a)  represent  the  set  of  all  adjacent  vertices  of  the 
vertex  represented  by  the  basic  feasible  solution  B  whose  cost  value  is  not  less  than 
that  of  B. 


To  Get  all  the  Basic  Feasible  Solutions  Representing  a  Degenerate  Vertex 

When  the  vertex  Vg  is  degenerate,  the  canonical  forms  of  all  the  basic  feasible  so¬ 
lutions  ,  which  represent  it,  may  be  obtained  by  iooking  at  the  canonical  form 

of  any  one  of  them  and  then  pivoting  among  the  non-zero  input-output  coefficients  in  the 
rows  corresponding  to  the  basic  variables  which  are  zero. 

Ranking  tne  Vertices  of  (l) 

Let  Sj,S2,...  be  a  ranking  of  the  basic  feasible  solutions  of  (l)  in  increasing 
order  of  Z  .  Suppose  we  already  know  the  basic  feasible  solutions  Sj ,S2> . . .  ^  in 

the  sequence.  It  is  intuitively  clear  that  the  next  element  in  the  sequence,  ,  must 
be  a  cost  nondecreasing  adjacent  vertex  of  one  of  *h a  vertices  represented  by  the  known 
basic  feasible  solutions  S| .... ,5^  j  .  We  shall  prove  this. 


Proposition  1:  Every  basic  feasible  solution  can  be  reached  by  taking  a  cost  nonde- 
creasing  path  from  through  the  vertices  of  (l)  . 


PROOF :  Consider  any  fasic  feasible  solu  ion  d  From  the  proof  of  the  simplex  algo- 
' i (hm  we  kn-^w  that  there  exists  a  cost  nonincreasing  path  moving  along  adjacent  vertices 
from  B  to  Sj  .  By  taking  the  same  path  in  the  reverse  direction  from  Sj  ,  we  reach 
B  from  Sj  by  moving  along  adjacent  vertices  along  a  cost  nondecreasing  path. 


Proposition  2:  Suppose  j  are  already  known.  Let  us  define 

p-1 

0  =  U  0{V<.  ]  -  fS. . Sn  ,)  p  =  2,3,... 

P  j_j  '  P-' 


(5) 
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Then  =  minimal  cost  solution  in  0^  . 

By  Proposition  I,  must  be  a  cost  nondecreasing  adjacent  vertex  of  one  of  the 

vertices  S| > . . . »S|c_j  .  But  is  the  minimal  cost  vertex  after  Sj,...,Sk  j  are 

excluded.  Hence  S.  =  minimal  cost  solution  in  0,  . 

k  k 

Now  the  algorithm  can  be  liven.  The  method  starts  with  the  finding  of  Sj  by 
the  simplex  method. 

General  Step:  Suppose  Sj . ^  have  already  been  obtained  and  we  are  trying  to 

find  out  Sk  .  Then  is  the  minimal  cost  basic  feasible  solution  among 

k-1  S.  S. 

U  {T.  1 :  j  such  that  x.^S.  and  C. 1  >  0}  -  fS. , . . . ,S.  . } 

;=i  J  J  1  J  l  k-l 


=  .U  0(S.)  -  fS, . sk-l J  .  (6) 

Of  course,  if  any  of  S.  are  degenerate  we  should  replace  0(S.)  by  0(V_  )  as  in 

'  '  5i 

equation  (4a). 

S. 

Thus  S.  can  be  easily  located  by  examining  the  values  Z (T . 1 )  for  i  *  l,...,k-l 
k  J 

S. 

and  j  such  that  x.^S.  and  C.  >  0  .  S,  is  that  new  basic  feasible  solution  in  (6) 
J  j  i  j  —  k 

which  is  distinct  from  Sj,...,Sk  j  and  which  has  least  cost  value  >  Zk  j  .  The  algl- 

rithm  is  stepwise  and  in  each  step  we  determine  an  additional  element  in  the  sequence 
of  ranked  vertices  Sj,  S^,... 


To  Organize  Computations:  Computationally,  this  may  be  done  by  storing  at: 

S. 

Array  1 :  All  the  Z(T.  )  values  for  all  S.  determined  so  far,  Vj  such  that 
Si  Si 

x.iS.  and  C.  >  0  and  T.  4  any  ' r  ‘he  known  S.  so  far  Of  course,  when  any  of 

J  '  J  —  J  1 

S.  are  degenerate,  we  should  scan  all  basic  feasible  solutions  which  represent  that 


same  vertex. 


Array  2:  All  the  basic  feasible  solutions  that  have  already  been  found  and  ranked. 
Each  of  the  S.'s  may  be  stored  in  terms  of  the  subscripts  of  the  basic  variables  in 
it,  arranged  in  increasing  order. 

S. 

Array  3:  The  basic  feasible  solutions  T.  corresponding  to  the  Z>values  stored 

in  Array  1.  Whenever  a  basic  feasible  solution  is  to  be  stroed,  store  the  subscripts  of 
the  basic  variables  in  it  in  increasing  order. 

It  is  convenient  to  locate  Array  I  and  Array  2  in  core  memory,  and  Array  3  on  tape. 

The  computations  required  to  get  the  next  element  in  the  sequence,  i.e.,  ,  are 

i)  to  scan  Array  1  completely  and  then  determine  the  least  value  there; 

ii)  to  identify  the  corresponding  basic  solution  from  Array  3-  This  is  . 

The  values  of  the  basic  variables  in  may  be  obtained  by  referring  to 

the  restrictions  (l).  If  it  is  required  to  find  out  some  more  elements  in 
the  sequence,  then 

iii)  delete  Z(Sk)  from  Array  1,  Sk  from  Array  3,  and  add  Sk  to  Array  2. 

iv)  find  out  the  canonical  form  of  and  all  its  cost  nondecreasing  adjacent 

vertices,  i.e.,  O(S^)  (or  0(V^  )  if  $k  is  degenerate).  Store  these 

ak 

basic  feasible  solutions  at  Array  3  and  their  Z-values  at  Array  1. 

If  the  problem  is  only  to  rank  all  basic  feasible  solution  for  which  Z  <  a  ,  then  we 
can  save  space  by  storing  in  Arrays  1  and  3  only  those  solutions  for  which  Z  <  a  . 

A  Numerical  Example:  We  apply  the  algorithm  to  the  following  fixed  cost  transportation 


2 


co»t  coefficient*  *  16. 


0,  -  00  -  67 


end  hence  It  may  be  nece**ary  to  rank  the  extreme  point*  of  the  trantportat ion  problem 
only  to  the  extent  of  l  <  22M*  ♦  67  ■  2281.  The  tableau  corre*pond I ng  to  2,  I*  given 
below. 


V* 

8  *‘J  ‘9 

2(T ,  V2270 


2230 

6  17 


2,  •  22  H*  .  0,-83 


r  •  2297 


Solution  S, 


Orly  Xjj-value*  correspond  I  ng  to  the  ba*ic  cell*  are  recorded  In  the  middle  of  the  cell. 

$. 

The  2(T.  )  value*  are  recorded  only  In  tho*e  r*onba»lc  cell*  where  it  >*  *  2281,  *lnce 


we  are  only  Interested  In  extreme  point*  which  have  2  <  2281. 

U*lng  equation  (6)  we  find  that  can  be  obtained  by  Introducing  x^ 


into  the 


{  .4 

CIJ 


18 

lit 

5 

21 

25 

6 

-2 

17 

7 

16 

8 

2k 

35 

6 

15 

31 

9 

40 

84 

29 

z2  •  2230 


02  •  <*6  f  •  2276 


Solution  S. 


Now  *  ®o  *  3®  «>vl  hence  it  i*  nece**ary  to  rank  the  extreme  point*  only  to  the  ex 
tent  that  Z  <  2230  ♦  30  •  2260.  In  the  tableau  corrctpond  i ng  to  S2  .  only  tho*e  non 

5, 

basic  cell*  which  lead  to  2(T/)  7260  have  been  marked. 

Using  equation  (6)  again,  we  find  that  Sj  I*  obtained  by  Introducing  lnto 


the  batl*  In  the  tableau  ot  S,  . 


J 


21 


It 


9  28 


•3 

6 

.  r 

|  * 

i  i6 

35 

6 

-3 

31 

3* 


79|  29 


REFERENCES 


r  I  ]  Hlrsch,  W.  H.  and  C.  6 .  Oantzig.  "Note*  on  Linear  Programming.  Part  XIX,  The 
Fixed  Charge  Problem,"  RAND  Research  Memorandum  No.  >38),  Santa  Monica,  Calif., 
095M. 

r2  ]  Hlrsch,  W.  N,  and  A.  J.  Hoffman,  "Extreme  Varieties,  Concave  Functions  and  the 
Flxad  Charge  Problem,"  Communlcat  Ions  on  Pure  and  Applied  Mathematics.  Vol.  XIV, 
No.  3  (August  1961). 

[3]  Bal Inski,  M.  1.,  "Fixed  Cost  Transports! Ion  Problems, "  Naval  Research  Logistics 
Quarterly.  Vol.  8.  No.  I  (March  I96I). 

[k]  Kuhn,  H.  W.  and  W.  J.  Baumol ,  “An  Approximate  Algorithm  for  the  Fixed  Charges 
Transports! Ion  Problem."  Naval  Research  Logistics  Quarterly,  (March  1982). 

[5]  Manna,  A.  S. .  "Plant  Location  under  Economies  of  Scale  Decentral Izat Ion  and  Com¬ 
putation  Management  Science.  Vol.  II,  No.  2  (November  198k). 


L 


